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ABSTRACT 


Observations of SN1006 have shown that ions and electrons in the plasma 
behind fast supernova remnant shock waves are far from equilibrium, with the 
electron temperature much lower than the proton temperature and ion tempera¬ 
tures approximately proportional to ion mass. In the ~360 km s“^ shock waves of 
the Cygnus Loop, on the other hand, electron and ion temperatures are roughly 
equal, and there is evidence that the oxygen kinetic temperature is not far from 
the proton temperature. In this paper we report observations of the He II A1640 
line and the C IV A1550 doublet in a 360 km s“^ shock in the Cygnus Loop. While 
the best £t kinetic temperatures are somewhat higher than the proton tempera¬ 
ture, the temperatures of He and C are consistent with the proton temperature 
and the upper limits are 0.5 and 0.3 times the mass-proportional temperatures, 
implying efficient thermal equilibration in this collisionless shock. The equilibra¬ 
tion of helium and hydrogen affects the conversion between proton temperatures 
determined from Ha line prohles and shock speeds, and that the efficient equili¬ 
bration found here reduces the shock speed estimates and the distance estimate 
to the Cygnus Loop of Medina et al. (2014) to about 800 pc. 

Subject headings: shock waves; ISM; supernova remnants; dust; ISM; individual 
(Cygnus Loop); ultraviolet; ISM instruments; HST(COS) 
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1. Introduction 


It is often assumed that astrophysical shock waves produce Maxwellian velocity 
distributions in the downstream plasma and that all particle species have the same 
temperature. However, in the modest number of cases where temperatures of different 
species can be measured, this is usually not the case. Shock waves in settings such as 
the solar wind, the interstellar medium, and galaxy clusters are collisionless, in that the 
shock transition occurs on a scale comparable to the proton gyroradius, which is orders of 
magnitude smaller than the collisonal mean free path. In this case, one might expect that 
a collisionless shock thermalizes a fraction of the kinetic energy of each incoming particle, 
leading to mass-proportional temperatures. The actual situation is more complex. 


Typical shock waves in the solar wind have modest Mach numbers, and they generally 
produce electron temperatures Te around 0.2 ti mes the proton temperatures Tp, while 


heavi er ions are preferentially heated, Tj > Tp (IGhavamian et ah 


20131: 


Korreck et ah 


20071). In the higher Mach number shocks in supernova remnants (SNRs), T^/Tp declines 
from about 1 to less than 0.1 as the shock speed increases from about 350 km s~^ to 


2000 km s ^ IGhavamian et ah 2001. 2002 

; 

dughes et ah 2000: Ghavamian et ah 2007a: 

van Adelsberg et ah 

2008; 

Medina et ah 

2014 

). Ion temperatures in SNR shocks have been 

more difficult to measure. 

but T,; ~ miTn/rrin for helium, carbon and oxvEen in a 2000-3000 


km s ^ shock in SN 1006 (iRavmond et ah 


Broersen et ah 


1995; 


Laming et ah 


1996 


Korreck et ah 


2004 : 


20131) . while the oxygen temperature is l ess than 1.7 times the proton 
temperature in a 350 km s“^ shock in the Cygnus Loop flRavmond et al.l 120031) . However, 


coulomb collisions affect the oxygen temperature in the Cygnus Loop. In the ~ 900 km s“^ 
shocks in the LMC supernova remnant DEM L71, the kinetic temperature of O derived from 
the O VI line widths is about 6.5 times the proton temperature, though this is somewhat 
uncertain due to interstellar absorption and the contribution of bulk motions to the line 
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width flGhavamian et alJl2007b|) . The shocks in galaxy clusters are as fast as the shocks 
in young SNRs, but the pre-shock temperatures are so high that the Mach number is only 
arou nd 3. Some seem to show rapid electron-ion temperature equilibration, while others do 


not flRussell et ah 


201211 ■ 


Direct m easurements of velocit ies in solar wind shocks generally show non-Maxwellian 


distributions fiThomsen et al. 


19901). Particle distributions in SNR shocks are non- 


Maxwellian in the sense that they have power-law tails of relativistic particles, but little 
is known about the velocity distribution of the vast majority of particles in the core of 
the velocity distribution. There is some evide nce for a non-thermal distribution at one 


position in Tycho’s SNR, but it is ambiguous flRavmond et al. 


201011 . Sir nulations with 


parti cle-in-cell codes tend to show Maxwellian cores with power-law tails fjCaprioli et al. 


2QM). 


In this paper we measure the kinetic temperatures of helium and carbon from prohles 
of the UV lines of He II and C IV obtained with the COS spectrograph on the Hubble Space 


Telescope. We observed a 360 


cm s ^ shock in th e Cygnus Loop for which electron-ion 


temperature ratio is close to 1 (iMedina et al. 


2014fl . In order to interpret the observations 


we needed the line prohle of the COS instrument with the G160M grating for an extended 
source, so we observed a position in the planetary nebula NGC 6853. We also compute the 
effects of optical depth and Coulomb collisions on the line prohles. We hnd that the helium 
and carbon kinetic temperatures, like the electron temperature, are close to the proton 
temperature. 
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2. Observations and Data Reduction 


The target position i n the Cygnus Loop co incides with the single orbit COS 


spectrum labeled Po s 1 in 


Medina et al 


Raymond et al. 


(120131) and the Ha observation labeled COSl in 


(1201411 ■ Nine more orbits were obtained with COS using the G160M grating 


set at 1577 A; 4 orbits on 2014 October 3 and 5 orbits on 2014 November 11. They were 
combined with the earlier exposure for a total observation time of 29,100 seconds. The 
circular entrance aperture is 2.5" in diameter. 

The po sition is shown i n Figu re [H which also shows the two positions behind the shock 


reported in 


Raymond et al. 


(I 2 OI 3 I) . It was chosen to be slightly behind the brightest Ha 


emission based on the ionization times required to reach He H and C IV in the shocked 
gas. The coordinates are RA(2000) = 20^* 54"* 43.61R, Dec(2000) = +32° 16' 3.53". For 
comparison, we include and X-ray image from Chandra and and NUV image from Galex. 
The X-rays show that the shock coincides with the jump in electron temperature. The 
structure in the Galex image is very similar to that seen in Ha. This is probably because 
the Galex NUV passband from about 1900 A to 2750 A is dominated by the 2-photon 
continuum of hydrogen in the this non-radiative shock, and that continuum from the 2s 
level of H is excited in a manner similar to the n=3 levels that produce Ha. 


The orientation of the dispersion axis of the COS instrument was close to the long 
direction of the hlament for the longer observations (18° and 16° away), though it was 
nearly perpendicular to the hlament (70°) for the earlier single orbit exposure. Because the 
dispersion direction was along the hlament for 90% of the exposure time, we assume that 
the emission hlled the aperture uniformly in the dispersion direction. 


m 


The position was cho sen because the huxes were kno wn from the s i ngle o rbit spectrum 


Raymond et al. 


(I 201 I and because the Ha prohle of 


Medina et al. 


(1201411 provides 


the proton temperature. It also provides an electron-proton temperature ratio of >0.8, 
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consistent with the determination of iGhavamian et al.l (120011) for a nearby position. In 
addition, the proper motions of the shocks in that re gion have been mea s ured and the 


electr on temperature is available from X-ray spectra flSalvesen et ah 


2009 


Katsuda et ah 


20081) . IR images from Spitzer and UV spectra from COS have been used to investigate 

tering of carbon atoms from grains behind this shock 


the destruction of dust and the spui 

derin 

(Sankrit et ahl 

2010 ; 

Ravmond et al. 

2OI3I 


20131) . Moreover, a FUSE spectrum of a position about 


2.7' to the NW al ong the same h 


O VI line prohle flRavmond et al 


amen t provides the oxygen kinetic temperature from the 


20031) . 


The expected line width is comparable to the instrumental prohle of COS for an 
extended source with the G160M grating, but no instrumen t al pro hle was available for 


France et al. 


( 2009 ) that the spectral 


extended sources other than the statement by 
resolution is about 200 km s“^. The prohle is not expected to be strictly Gaussian, but to 
be ahected by the shape of the entrance aperture. To obtain the instrument prohle for an 
extended source of narrow emission lines through the COS aperture, we observed a position 


in the p 


anetary nebula NCG 6853 (The Dumbbell) 28" S and 2" E of the central star. 


Goudis et al. 


Barken (11984 ) repo rted emission line huxes from an lUE spectrum at that position, and 


( 1978 ) show [O III] line prohles of NGC6853. The two positions closest to the 


GOS position have line widths (FWHM) of about 30 and 60 km s“^. An archival STIS 
spectrum of the central star of NGG 6853 shows saturated absorption against the stellar 
continuum in the G IV lines between 0 and -40 km s“^, which probably removes the blue 
wings of the emission lines, so the width is about 30 km s“^, which is small compared to 
the GOS instrumental FWHM. At worst, the FWHM of the emission lines could be about 
40 km s“^, which would mean that we overestimate the width of the instrumental prohle by 
4%. We hnd that the prohles of the G IV A1550 doublet and the He H A1640 line in NGG 
6853 can be ht with pairs of Gaussians. For the G IV lines, the FWHM is 0.40 A and the 
separation is 0.86 A, while for the He H line the best ht is FWHM = 0.47 A and separation 










































7 


= 0.93 A. The larger width of the He II line is at least partly due to the fact that He H 
A1640 is an unresolved multiplet with components spread over about 0.16 A. 


Figure [2] shows the observed line profiles along with the instrument profiles from the 
planetary nebula measurements described above. The Cygnus Loop profiles are somewhat 
wider than the instrument profile, but not by much. We fit the profiles with Gaussians 
convolved with the instrument profiles described above. Figure [3] displays these hts, and 
Table 1 shows the results with 1-a errors, along with the corresp ondi ng kinetic temperature s 


The proton and O VI temperatures from 


Medina et ah 


(120141 ) and 


R aymond et ah 


(1200311 


are included for comparison. Note that the proton temperature is not obtained directly 
from the Ha line width. Charge transfer between protons and neutral atoms produces 
a population of neutrals given by the proton thermal d istribution convolved with t he 


charge transfer cross section times the relative velocity (j Chevalier fc R,aymondlll978l ). so 


the population of fast neutrals is similar to, but not th e same as, the proton velocity 


distribution. We use the model of 


Chevalier et ahl (1198011 to derive the proton temperature. 


3. Analysis 

Table 1 shows that the kinetic temperatures of helium and carbon are indistinguishable 
from the proton temperature, with upper limits of 1.4 and 2.9 times the proton temperature, 
or 0.35 and 0.24 times the mass-proportional temperatures, respectively. Here we discuss 
corrections for optical depth and Coulomb scattering used to obtain those limits. 


3.1. Optical depth effects 


The optical hlaments of the Cygnus Loop and other SNRs are t angencies of the 


line of sight to a thin, rippled sheet of emitting plasma ((Hester 


19871) . Thus they are 

















essentially thin slabs of emission seen edge-on. Resonance lines snch as the C IV donblet 
can have snbstantial optical depths along the line of sight, bnt small optical d epths in the 


perpendicular direct 

ion, s 

1992; 

Cornett et al. 

1992) 


don, s o that photons are scattered ont of the line of sight flLong et al 


19921) . The net effect for lines of sight near tangency is to rednce the 


i ntensity and, because the optical depth is highest at line center, to broaden the profile 


(iRavmond et al. 


299311 . 


The ratios of the intrinsic emissivities and of the scattering cross sections for the C 
IV doublet are both 2:1, but the observed intensity ratio is 1.65:1. In the single scattering 
limit, the intensity at any wavelength is proportional to (l-e“’’), and the observed intensity 
ratio corresponds to optical depths at line center of about 9.65 and 9.33 for A1548 and 
A1559, respectively. Those optical depths imply increases of the FWHM by 12% and 5.3%, 
respectively. We correct the carbon kinetic temperatures for these broadenings in Table 1. 
We note that the A1548 line is nominally wider than the A1559 line, though the difference 
is easily within the uncertainties. If the apparent larger width of A1548 were attributed 
to scattering, the observed relative intensities would be close to 1:1. The optical depth 
correction for the width of the O VI line is more sev ere because of its larger optical depth. 


and we use the value given by 


Raymond et al. 


(1299311 . 


The He II A1649 line is optically thin, and it is not directly affected. However, some 
He H A256 photons are absorbed and 12% of the absorbed A256 photons are converted into 
A1649 photons. Since the optical depth is largest near line center, more absorptions occur 
close to line center, producing excess A1649 photons near line center and making the line 
narrower. However, the A256 optical depth perpendicular to the shock is about 6, and much 
larger along the line of sight. Therefore, He H Ly/1 photons are converted to A1649 photons 
over much of the line profile, and we will ignore this effect on the line width. 
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3.2. Coulomb Equilibration 


We are interested in the ion kinetic temperatnres jnst behind the shock front, bnt it 
takes some time for carbon to be ionized to C IV, and both He II and C IV emit over a 
range of distances behind the shock determined by the ionization times. Dnring those times. 
Coulomb collisions will bring the kinetic temperatures toward equilibrium with the proton 
temperature. Thus we must compute the changes in the ion temperatures due to Coulomb 
equilibration to see how accurately the observed kinetic temperatures reflect the values at 
the shock. Both because the ionization rates are smaller for higher ionization states and 
because the Coulomb collision rate scales as we expect that this equilibration will be 
negligible for He H, larger for C IV and quite signihcant for O VI. 


We compute the time-dependent ionization states of He, C and O using the electron 
temperature from X-rays of 2x10^ K with t he ionization rates of iDerel (120071) as compiled 


in version 7 of CHIANTI ( Landi et ah 


20121) . We use a prot on temperatur e of 1.8 x 10® K 


from the Ha prohle with ion-proton equilibration rates from 


Spitzerl (119561) . The difference 


between the two assumed temperatures lies within the uncertainties. As a limiting case, we 
assume mass-proportional initial temperatures, Ths = 7.2 x 10® K, Tc = 2.2 x 10^ K and 
To = 2.9 X 10^ K, with all the ions starting in the singly ionized state. The density does not 
matter, because the ionization and Coulomb collision rates scale in the same way with Ug- 


The top panel of Figure 0] shows the kinetic temperatures of He H, C IV and O 
VI ions as functions of time for a density of 2.0 cm“^. It is apparent that carbon and 
oxygen approach equilibrium relatively quickly because of the dependence of the 
Coulomb collision rates. This plot is not very useful, however, because nearly all of the 
oxygen is in lower ionization states at short times and in higher ones at long times. The 
middle panel shows the ionization fractions of He H, C IV and O VI as functions of the 
kinetic temperatures of those ions. The result, as expected, is that the average kinetic 
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temperature of the He II is the same as the initial temperature. The carbon ionization 
fraction peaks when the C IV kinetic temperature had dropped noticeably, and the O VI 
kinetic temperature has dropped even more sharply before the peak ionization fraction has 
been reached. 


The bottom panel in Figure 0] shows the fractional contribution to the emission for 
each ion as a function of temperature in 10® K bins. The emission of He H is conhned to a 
single bin, while the emission of C IV is concentrated between 1.6 x 10^ and 2.1 x 10^ K, 
with a tail at lower temperatures. The O VI emission peaks in the bin between 2 x 10® and 
3 X 10® K with a long tail at higher temperatures. The peak of the O VI contribution in 
the low temperature bin close to the proton temperature comes about because the oxygen 
spends a lot of time at those kinetic temperatures, which compensates for the modest 
ionization fraction. The average temperatures are 7.2 x 10®, 1.7 x 10^ and 5.6 x 10® K 
for He H, C IV and O VI, respectively, or 1.0, 0.77 and 0.2 times the mass-proportional 
temperatures. We computed similar models for a range of intial temperatures and hnd that 
the upper limits in Table 1 are compatible with intial temperatures of 0.5, 0.3 and 0.5 times 
the mass proportional temperatures for He, C and O, respectively. Smaller values would be 
needed to match the best fit temperatures. 


The fact that C IV and O VI are formed over a range of temperature means that their 
line prohles are the sum of Maxwellians, so they will deviate from Gaussians. The C IV line 
is formed over a range of about 25% in temperature or 13% in thermal speed, so detecting 
the departure from Maxwellian would require better data than the spectrum presented 
here. The O VI profile from a shock with mass-proportional heating would give a core with 


a temperature just above the proton temperature p 
of the flux. While the some of the prohles shown in 


us broad wings 


Raymond et al. 


contai ning a large part 


(120031 ) show a hint of a 


wing on the red side, there is no indication of the strong wings predicted by a model with 
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mass-proportional temperatures. 


Thus we can ignore Coulomb collisions for the He II line, a nd the C IV lines wid 
are affected at the 12% level, but the O VI prohles presented by 


severely affected . In the faster shocks in D EM L71 (IGhavamian et ah 


(iRavmond et ah 


1995 


Korreck et ah 


Raymond et ah 


fl2003l) are 


ihs 


2nn7bh and SN1006 


2004j) . the higher electron and ion temperatures 


greatly reduce the ionization time scales and increase the Coulomb collision time scale, so 
that Coulomb collisions can be neglected even for O VI. 


3.3. Interpretation 


First we discuss the kinetic temperatures that might be expected in a collisionless 
shock. To hrst order, such a shock simply thermalizes 3/4 of the shock speed independently 
for each species, leading to mass- proportional ion temperatures, approximat ely as seen in 


the 3000 km s ^shock in SN1006 fiRaymond et al 


1995; 


Korreck et al. 


200411 . 


The situation is more complicated for species that enter the shock as neutral atoms 
or in dust grains. Mu ch of the helium is ne utral, and about 75% of the carbon is in the 


form of PAHs or dust fiRaymond et al. 


20131) . The PAHs dissociate very quickly behind the 


shock and contribute to the observed line prohle flMicelotta et al 


2010 l) . while dust grains 


sputter away more gr adually over tens of arcseconds at the distance of the Cygnus Loop 


fiRaymond et al. 


20131) and do not contribute much to the prohle at the shock. Neutrals that 


pass through the shock and become ionized downstream are expected to behave as pickup 


ions fiRaymond et al. 


200811 . For quasi-perpendicular shocks, the velocity component along 


the held is conserved while the component perpendicular to the held forms a ring beam 
in velocity space. The unstable ring distribution quickly scatters into a hollow bispheric al 


shell in velocity space by emitting and absorbing Alfven waves fjWilliams fc Zank 


19941: 
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Isenberdll995l) . The distribution then scatters into a Maxwellian, though the temperatures 
parallel and perpendicular to the held may differ. If the distribution is too anisotropic, 
either firehose or mi rror instabilities will drive the anisotropy toward marginal instability 


(iMaruca et ah 


20121) . Such a pickup ion process in the partially neutral hydrogen provides 


one of several possible interpretatio ns for the non-Iv 


’g” in Tycho’s supernova remnant (iRavmond et ah 


axwe 


lian Ha profile observed at knot 


20101 ) ■ and it is responsible for the 


motion of oxygen ions along field lines seen in images of the sungrazing comet C/2011 W3 


(Lovejoy) flRavmond et af 


20141) . The poor fit of the He H profile to a Gaussian velocity 


distribution, = 2.4, might suggest a pickup ion contribution added to a Maxwellian, but 
the data do not warrant a strong conclusion. We note that the sharp peak to the right of 
the centroid is much narrower than the instrumental resolution and must therefore be a 
result of the modest number of counts in the line. 

The comet observations also confirm that a significant fraction of the energy in the 
ring distribution can be lost to waves, depending on th e Alfven speed and angle between 


the flow and the magnetic field (j Williams fc Zank 


19941 ). For the shock observed here, the 


post-shock Alfven spe ed is around 20-30 km s and the loss to waves should be small 


(IRaymond et ah 


20081) . The final temperature of the pickup ions depends on the angle 


between the field and the shock, and for quasi-perpendicular shocks it will be between 0.5 
and 1 times the mass-proportional temperature. 

Many of the C and O atoms are already ionized when they pass through the shock. 
Because of their high mass to charge ratios, these ions are not slowe d as much by the 


electro static field in the shock wave as are the protons. According to 


( 1997 ) they will be about 25% hotter than the protons. Ihee &: Wul ( 2000 ) and 


Fuselier fc Schmidt 


Zimbardo 


(120111 ) also present models for strong preferential heating of heavy ions in shocks. Those 
models predict very high ratios of perpendicular to parallel temperatures, T^/T\\. Depending 
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upon the angle between the magnetic held and our line of sight, we could in principle 
observe mainly Ty and therefore a low temperature. However, large 7±/7"i| would generate 
mirror mod e instabilities and, for the high plasma /3 in this shock, sharply reduce the 


anisotropy (iMaruca et ah 


201211 ■ Solar wind observations of modest Mach number shocks 


show some preferential heating of heavy ions, but there is little correlation with shock 


parameters flKorreck et ah 


2003 ). 


There is also transfer of energy among species by plasma waves and turbulence at 
the shock. Waves can be excited upstream of the shock jump by str eaming instabilities 


generated eit 


ler by particles rehect ed by the shock or by cos mic r avs flThomsen et a1 


Laming et al 


201411 . For example. 


Ghavamian et al. 


(2007a) and 


Rakowski et al. 


1990 


fcoosh 


show that lower hybrid waves in a cosmic ray precursor can heat the electrons to T^. ~ Tp in 
shocks at about the speed of the one observed here. 


Our observations give up per limits to He anc 
proton temperature found by 


Medina et al 


C temperatures of 1.4 and 2.9 times the 


(12014^ . After considering Coulomb collisions. 


the upper limits on the initial tempera 


and the O VI profile of 


Raymond et al. 


ures a re 2 and 4 times the proton temperature. 


(1200311 gives an upper limit of 8 times the proton 


temperature. This contrasts with the nearly mass proportional values (4, 12 and 16) seen 

in the 2000-3000 km s ^.shock in SN1006 (Raymond et al. 

1995: 

Korreck et al. 2004). and 

it is similar to the indications of nearly complete electron-ion ec 

uilibration in the Cvenus 

Loop compared with T^/T, ~ 0.05 in SN1006 (Ghavamian et ah 

2002 : van Adelsberg et al. 


20081) . It indicates rapid temperature equilibration by plasma waves and turbulence in the 
Mach 35 Cygnus Loop shock, while the equilibration in the Mach 200-300 SN1006 shock is 
ineffective. 

The effective equilibration of the particle species has an i mportant implication fo r the 


derivation of shock speeds from Ha line widths. As shown by 


Chevalier et al 


(1198011 and 
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van Adelsberg et al. 


(120081) . if the energy dissipated in the shock is shared among electrons 


and ions, a higher shock speed is needed to account for a given line width. On the other 
hand, if helium is e fficiently equilibrated , the proton temperature behind a shock of a given 


speed is increased. 


Medina et al. 


(120141) assumed that helium remained mostly neutral 


through the region where the Ha line is formed and that the Tg and Tj are equal, so that 
the post-shock temperature is given by 1/2 the thermalized proton speed. On the other 
hand, we hnd that helium is brought to close to the ion temperature close to the shock, in 


which case the kinetic energy of the helium is shared with the protons. While 


Medina et ah 


(120141) found a shock speed of 405 km s“^, a shock speed as low as ~ 360 km s“^could 
account for the Ha prohle if Th = The where the Ha line forms. (We assume a helium 
abundance of 10% by number. If most of the helium is neutral and it is mostly ionized only 
after the hydrogen is ionized, the effect will be smaller.) When that lower shock speed is 


combined with the proper motion measured by 


Salvesen et al. 


( 20091) . the distance to the 


Cygnus Loop is reduced from 890 pc to 800 pc. That partially alleviates the disagreement 
between the proper motion distanc e and the upper limit of 640 pc to the distance from the 


detection of SNR absorption lines (iBlair et al 


20091) in the spectrum of a background star. 


but it still leaves a substantial discrepancy. Of course, if a signihcant fraction of the shock 
energy goes into cosm ic rays, a correspondingly higher shock speed would be required. 


Shimoda et al. 


(120151 ) point out that comparison of the shock speed derived from the 


proper motion with the shock parameters determined from the shock jump conditions can 


be misleading, as a shock in a realistic ISM density di stribution is s o mew 
most positions. However, for the shock observed here. 


Medina et al. 


lat oblique at 


(12014^ measured the 

offset between the velocity centroids of the broad and narrow Ha to be only 9 km s“^, so 
the obliquity is very small. 
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4. Summary 

We have measured He and C kinetic temperatures behind a moderate velocity shock 
in the Cygnus Loop. Unlike the fast shock in SN1006, where the temperatures are 
mass-proportional, in the Cygnus Loop they are nearly the same. We have shown that 
Coulomb equilibration is not important for the He H line, and it affects the C IV lines 
at the 33% level, so the equilibration must occur at the shock front by means of plasma 
turbulence. However, Coulomb collisions severely affect the kinetic temperature of O VI if 
it is initially higher than the proton temperature, so that the C IV and He H lines provide 
better constraints on temperatures at the shock than do the O VI lines for this relatively 
slow SNR shock. The rapid equilibration of helium suggests that the post-shock protons 
are partly heated by helium, which would reduce the shock speed required, and hence the 
distance inferred from proper motions, by as much as 11%. It is important that the the 
thermal equilibration of protons, electrons and helium be considered when converting Ha 
line widths of Balmer line hlaments into shock speeds. 

The degree of ion-ion equilibration and the related heating processes may be important 
for determining the efficiency of injection of the different species into the cosmic ray 
acceleration process. Besides that, it is potentially useful as a diagnostic for the plasma 
waves that produce the shock transition in a collisionless shock. 

We determined the instrumental prohle of the COS spectrograph for an extended 
source with the G160M grating by observing a planetary nebula, and this may be useful 
for other studies of extended sources that £11 the COS aperture. The ~ 200 km s“^ 
wide instrumental profile can be described as a pair of Gaussians with FWHM = 0.40 A 
separated by 0.86 A. 


The authors thank the referee for requesting an additional plot, which revealed a 



numerical error in the average temperatures, and for other useful suggestions. This work 
was performed under grants HST-GO-12885 and HST-GO-13436 to the Smithsonian 
Astrophysical Observatory. RJE was supported by NASA contract NAS8-03060. JR and 
PG thank Lorentz Genter Workshop on Particle Acceleration from the Solar System to 
Gosmology for useful discussions. 
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Table 1 

Fit Parameters 

Temperatures in units of 10® K 


Line 

FWHM (km s"^) 

Tkin (10® K) 

T,orr (10® K) 
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2541}®“ 
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1 + 
o o 

CO CO 



He H A1640 


2.011® 

2.011® 

2.46 

C IV A1548 

C IV A1550 

106125 

^79+35 

' ^-71 

K Q+3-2 

o.y_3 7 

9 7+3.2 
z,. 1 _2 7 

4 - 711 ^ 

+-n_2,4 

1.33 

1.43 

0 VI A1032 

120 

5.R 

<3.2^= 



a Ha line width from 


Medina et ah 


(12014 


b Hydrogen temperat ure based on mode 


after ch arge transfer ((Medina et al. 


c From 


Raymond et ah 


2oi4 


proton temperature leading to observed line width 


(120031) . The temperature is 3.2 x 10® K after correction for optical 


depth effects. 





















21 



Fig. 1.— COS aperture positions of 


Raymond et al. 


(120131 ) overlaid on Ha, Chandra X-ray 


and Galex NUV images. The lower right hand panel shows the apertures overlaid on a 3 
color superposition of Ha (red), NUV (blue) and X-rays (green). The position observed here 
is the one at the upper left on the Ha filament. 
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Fig. 2.— Observed C IV doublet and He II A1640 lines from the Cygnus Loop (solid) 
compared with the instrument prohles derived from the spectrum of the planetary nebula 
NGC 6853 (dashed). 
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Fig. 3.— Observed C IV doublet and He H A1640 observed profiles (solid) compared with 
best £t Gaussians convolved with the planetary nebula line profiles (dashed). 
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Fig. 4.— Top panel. Kinetic temperatures of He II, C IV and O VI as functions of time 
behind the shock front. Middle panel. Ionization fractions of He H, C IV and O VI as 
functions of kinetic temperature. While Coulomb collisions have no effect on helium before 
He H is ionized away, and only a small effect on the the temperature of carbon before C IV 
is ionized to C V, the temperature of oxygen is driven toward the proton temperature by 
the time O VI is ionized. Bottom panel. Contribution fractions to the emission of each line 
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